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(Dated: February 7, 2008) 

We present a measurement of the cross section of the process e + e~ — > tp(2S) from threshold 

up to 8 GeV center-of-mass energy using events containing initial-state radiation, produced at the 
PEP-II e + e" storage rings. The study is based on 298 fb 1 of data recorded with the BABAR 
detector. A structure is observed in the cross-section not far above threshold, near 4.32 GeV. We 
also investigate the compatibility of this structure with the F(4260) previously reported by this 
experiment. 

PACS numbers: 14.40.Gx, 13.25. Gv, 13.66. Be 



Until recently, charmonium spectroscopy has been well 
described by potential models. Observations of the 
X(3872) [H and the F(4260) decaying into ir+ij- J/t/j 
complicate this picture, and have stimulated both ex- 
perimental and theoretical interest in this area. The 
y(4260) can be produced by direct e + e~ annihilation 
and is therefore known to have J PC = Weak 
evidence for the Y (4260) structure in B decays was 
also reported by BABAR [3]. In addition, the F(4260) 
has been confirmed by the CLEO-c experiment in di- 
rect e+e- -> Y"(4260) interactions where the F(4260) 
is detected in decays to ir + ir~ J/tp and 7r°7r°J/^ |4; the 
observation of the latter mode and the measured ratio 
£(y(4260)^ ir°ir°J/tp)/B{Y (4260)-+ n+ir-J/ip) « 0.5 
implies that the F(4260) has isospin zero, as expected 
for a charmonium state. 

It is peculiar that the F(4260) is wide and yet has 
a large branching fraction into the hidden charm mode 
7r + 7r~ J/tp , and that at the Y (4260) mass the cross section 
for e + e~ — > hadrons exhibits a local minimum. Many 
theoretical interpretations for the F(4260) have been 
proposed, including unconventional scenarios: quark- 
antiquark gluon hybrids || and hadronic molecules @. 
We undertook this study with the intent of clarifying the 
nature of the F(4260). 

In this Letter we study the process e + e~ — > 
■K + 7r~tp(2S), tp(2S) — > 7r + n~J/tp, for e + e~ center-of- 
mass (CM) energies from threshold up to 8 GeV using 
initial-state radiation (ISR) events. The ISR cross sec- 
tion for a particular hadronic final state / is given by 

*2^ = W(s,x).*Ml-x)), (1) 

where s is the square of the e + e~ CM energy, x = 
2E 1 /y/s is the ratio of the photon energy to the beam 
energy in the e + e~ CM frame, and W(s,x) is the spec- 
trum for ISR photon emission for which we use a calcu- 
lation good to 0(a 2 ); the effective CM energy \fp is the 
invariant mass of the final state m = ys(l — x). 

We use data recorded with the BABAR detector [7J at 
the PEP-II asymmetric-energy e + e~ storage rings, lo- 
cated at the Stanford Linear Accelerator Center. These 
data represent an integrated luminosity of 272 fb -1 
recorded at y/s — 10.58 GeV, near the T(4S) resonance, 
and 26 fb -1 recorded near 10.54 GeV. 

Charged-particle momenta are measured in a tracking 



system consisting of a five-layer double-sided silicon ver- 
tex tracker (SVT) and a 40-layer central drift chamber 
(DCH), both situated in a 1.5-T axial magnetic field. 
An internally reflecting ring-imaging Cherenkov detec- 
tor (DIRC) provides charged-particle identification. A 
CsI(Tl) electromagnetic calorimeter (EMC) is used to 
detect and identify photons and electrons, while muons 
are identified in the instrumented magnetic-flux return 
system (IFR). 

Optimized selection criteria are chosen based on a sim- 
ulated sample of e + e~ — > 7 ISR 7^ + 7^ - '0(2«S , ) events and a 
sample of e + e~ — > j ISB tp(2S), tp(2S) — > 7r + 7r~ J/tp candi- 
dates in data, which serves as a clean control sample [f| . 

A candidate J/tp meson is reconstructed via its decay 
to e + e~ or fi + fj,~ . The lepton tracks must be well re- 
constructed, and at least one must be identified as an 
electron or a muon candidate. An algorithm to recover 
energy lost to bremsstrahlung is applied to electron can- 
didates. An e + e~ pair with its invariant mass within the 
interval of (-100,-1-40) MeV/c 2 of the nominal J/tp mass is 
taken as a J/tp candidate. For a /x + fj,~ pair, the interval is 
(-60,-1-40) MeV/c 2 . The J/tp candidate is then kinemati- 
cally constrained to the nominal J/tp mass and combined 
with a pair of oppositely-charged tracks identified as pion 
candidates. The 7r + 7r~ J/tp combinations with invariant 
mass within 10 MeV/c 2 of the nominal tp(2S) mass are 
taken as tp{2S) candidates. Another pair of oppositely- 
charged pion candidates (primary pions) is then com- 
bined with the tp{2S) candidate. The 7r + 7r _ V'(2S') mass- 
resolution function is well described by a Cauchy distri- 
bution [9] with a FWHM of about 7 MeV/c 2 . We do 
not require observation of the ISR photon (7 LSK ) as it is 
preferentially produced along the beam directions. 

We select e + e~ — > 7 ISR 7r + 7r _ -0(2iS') events with the fol- 
lowing criteria: (1) there must be no additional well- 
reconstructed charged tracks in the event; (2) there must 
be no well-reconstructed 7r° or rj — > 77 in the event; (3) 
the transverse component of the visible momentum in the 
e + e~ CM frame, including that of the 7 ISK when it is re- 
constructed, must be less than 2 GeV/c; (4) the difference 
(Ap*) between the measured it + i{~tp(2S) momentum 
and the value expected for it in an ISR 7r + 7r~ tp(2S) event, 
that is, (s — m 2 )/(2y/s), where m is the n + n~tp(2S) in- 
variant mass, must be within [—0.10, +0.06] GeV/c; (5) 
cos 81, where 9i is the angle between the lepton i + mo- 
mentum in the J/tp rest frame and the J/tp momentum 
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in the e + e~ CM frame, must satisfy | cos^| < 0.90; and 
(6) the invariant mass of the ir + ir~ pair in ip(2S) decay 
must be greater than 0.4 GeV/c 2 in order to suppress the 
combinatorial 7r + 7r _ J/0 background. 

A clean ip{2S) signal is apparent in Fig. [TJ An ex- 
amination of the ir + Tr~ip(2S) combinations reveals that 
about half the background results from recombinations 
within the same 2(ir + Tr~)J/ip system where at least 
one of the primary pions is combined with the J/tp 
to form a tt + tt~J/iP candidate. After subtracting the 
self-combinatorial background, we estimate 3.8±1.1 non- 
ij)(2S) background events in the final sample of 78 events 
within the 0(25) mass window. 
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FIG. 1: The invariant mass distribution for all ir + n~J/ip 
candidates where more than one entry per event is allowed. 
The solid curve is a fit to the distribution in which the ij}(2S) 
signal is described by a Cauchy function and the background 
by a quadratic function (represented by the dashed curve). 
The arrows indicate the ip(2S) mass window. 

In Fig.[5]the distributions of (a) Ap* and (b) cos 9* for 
2(7r + 7r~) J/ip candidates, where 8* is the angle between 
the positron beam and the (ir + Tr~ir + TT~ J/0) momentum 
in the e + e~ CM frame, are shown and compared to ex- 
pectations from simulations. There are 16 events that 
have a well reconstructed gamma with energy greater 
than 3 GeV, while the Monte Carlo simulation predicts 
16.4 for the same total number of ISR 7r + 7r~0(25) can- 
didates. Furthermore, all events within |cos6**| < 0.9 
are accompanied by a reconstructed gamma with energy 
greater than 3.0 GeV. We find excellent agreement in the 
ISR characteristics between the data and signal Monte 
Carlo sample. The good agreement in the Ap* distri- 
bution rules out any significant feeddown from higher 
mass charmonia decaying to the tp{2S) with one or more 
undetected particles. As an example, the Ap* distri- 
bution for -0(4415) — > ir + ir~Tr°ip(2S) events would peak 
around —0.2 GeV/c with a long tail extending to well be- 
low -0.2GcV/c. We estimate the non-ISR 7r+7r^0(2S") 
background to be less than 1 event. 

The track quality, particle identification information, 
and kinematic variables of all pion candidates are exam- 
ined, and displays of the events are scanned visually to 
check for possible track duplications and other potential 
problems. No evidence for improper reconstruction or 




0.05 
Ap* (GeV/c) 



FIG. 2: The distributions of (a) Ap* and (b) cos#* of the 
2(ir + iv~)J/ip combination in the e + e~ CM frame are shown 
for data (solid dots) and Monte Carlo simulation of the signal 
(histogram) normalized to the total number of the observed 
data events. 



event quality problems is found. 

The 2(tt + 7t~)J/iP invariant-mass spectrum up to 
5.7 GeV/c 2 for the final sample is represented as data 
points in Fig. [3] A structure around 4.32 GeV/c 2 is ob- 
served in the mass spectrum. 

To clarify the peaking structure observed in Fig. [3J we 
perform an unbinned maximum likelihood fit to the mass 
spectrum up to 5.7 GeV/c 2 in terms of a single resonance 
with the following probability density function (PDF): 



P(m) =Na- e(m) ■ (w(s, x) ■ 2m/ s] £ 

M 2 ■ r M • I> • ($(m)/$(M)) 



(2) 



( M 2- TO 2)2 + ( Mrtot) S 



+ B{m) , 



where M, T tot , T ee ,Tt,N are the nominal mass, to- 
tal width, partial width to e + e~, partial width to 
7r + 7r _ i/;(2S'), and yield for a resonance, respectively, and 
m is the 2(7r + 7r~) J/0 invariant mass, e(m) is the mass- 
dependent efficiency, $(m) is the mass-dependent phase- 
space factor for a S'-wave three-body ir + ir~ip(2S) sys- 
tem, a is a normalization factor, and B(m) is the PDF 
(the shaded histogram in Fig. [3]) for the non-0025 1 ) back- 
ground. The shape of B was obtained from ip(2S) side- 
band events with its integral fixed to 3.1 events corre- 
sponding to the mass region in the fit, where the total 
number of events is 68. The mass dependence of r tot is 
ignored in the fit. 

We perform fits to the distribution in Fig. [3] to test 
hypotheses that the data are a result of the decay of 
the T(4260) (dashed curve) using resonance parameters 
fixed to those of Ref , and alternatively those of the 
-0(4415) (not shown) with the mass and width taken from 
Ref [10(. In the third fit (solid curve) we assume a sin- 
gle resonance whose mass and width are free parame- 
ters, which are then found to be (4324 ±24) MeV/c 2 and 
(172 ± 33) MeV (after unfolding mass-resolution) by the 
fit. We calculate the % 2 /dof value for each fit to test these 
hypotheses. In the calculation the events in Fig. [3] are re- 
grouped so that at least seven events are expected in each 
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FIG. 3: The 2(tv + tv )J/ip invariant mass spectrum up to region m(2(7r + 7r") J/ip ) < 5.7 GeV/c 2 for the final sample. 

5.7 GeV/c 2 for the final sample. The shaded histogram rep- Only one entry per event is included in the plot, as described 

resents the fixed background and the curves represent the fits m the text. The histogram shows the distribution for Monte 

to the data (see text). Carlo events (see text). 



bin. The x 2 /dof values are found to be 21.3/8, 54.4/7, 
and 7.3/7 for hypotheses of the F(4260), the ^(4415), 
and a new resonance, respectively, corresponding to x 2 - 
probabilities of 6.5 x 1CT 3 , 2.0 x 1CT 9 , and 29%. The 
low probabilities associated with the F(4260) and the 
ip(4A15) indicate that the structure is not consistent with 
the ?/>(4415), and is not well described by the y(4260) ei- 
ther. We also perform a fit including both the y(4260) 
and ■0(4415) plus their interference, and find the x 2 /dof 
value to be 17.8/6, corresponding to a x 2 -probability of 
6.7 x 10~ 3 , but no much improvement from the fit to the 
Y"(4260) only. In order to further compare the structure 
reported here with the F(4260) reported in Ref. 0], we 
perform simultaneous fits to both the ir + Tr~ip(2S) mass 
spectrum in Fig.[3]and the tt + tt~ J/tp mass distribution in 
[2J under the hypotheses that (1) both signals are a single 
resonance and (2) these signals are manifestations of two 
independent resonances, with a single resonance for each 
signal. The PDF as used in Ref. is applied to the fit to 
the 7r + 7r - J/ip mass distribution. The logarithmic likeli- 
hood obtained from the single-resonance hypothesis (1) is 

5.4 units less than that obtained from the two-resonance 
hypothesis (2), which corresponds to a x 2 -probability of 

4.5 x 1CP 3 for the single-resonance hypothesis assuming a 
X 2 distribution for the difference in the logarithmic likeli- 
hood between the two hypotheses. However, none of the 
probabilities associated with the Y"(4260) can exclude the 
possibility that the structure observed is a manifestation 
of a new decay mode for the F(4260). 

The primary tt + tt~ invariant mass distribution for the 
selected events within m(2(7r + 7r~) J/0 ) < 5.7 GeV/c 2 is 
shown in Fig. |4j For the two events having more than 
one ip(2S) candidates, the dipion invariant mass is only 
included for the ip(2S) candidate closest to its nominal 
mass. The Monte Carlo distribution is also shown in 
Fig-Hlfor a single resonance decaying to it + it~iI)(2S) in a 
S-wave three-body phase-space using the resonance pa- 



We extract the energy-dependent cross section for 
e + e~ — > 7r + 7r~i/;(2S) up to 8 GeV for the final sample. 
The average cross section over a mass range of width Am 
is calculated as 

m+Am/2 



a(m) = J <J ( a; ) dx j Am 

J m— Am/2 

1 x 



C-B- Am^\2mi/s- W(s,l-m%/s) ■ 



(3) 



where C is the integrated luminosity, B is the product 
of B(ip(2S) -* tt+tt-J/^) and B(J/i/j -* £+T), the 
sum is over all events within the mass range, rrii is the 
2(7r + 7r _ ) J/ip invariant mass, and S{ is the estimated effi- 
ciency at that mass. The measured cross section is shown 



in Fig. [5] and the numerical results can be found in llj , 
where the background has been subtracted from bins 
with non-zero content. The energy-dependent selection 
efficiency (solid histogram in Fig. [5]) is determined from 
Monte Carlo events for which the ip(2S) polarization 
has been properly considered while the primary ir + ir~ 
is generated in S'-wave phase-space. The uncertainty in 
the selection efficiency due to model dependence is es- 
timated from the efficiency difference between S'-wave 
phase-space model and multipole model [12| in the pri- 



mary 7r + 7r _ generation. The main systematic uncertain- 
ties are listed in Table HI and are added in quadrature, 
resulting in a total systematic uncertainty of 12.3%. 

In summary, we have used ISR events to study the ex- 
clusive process e + e~ — * 7r + 7r~-0(2S) and to measure its 
energy-dependent cross section from threshold to 8 GeV 
CM energy. A structure is observed at ~ 4.32 GeV/c 2 in 
the 7r + 7r~'0(2S) invariant mass spectrum that is not con- 
sistent with the decay ?/>(4415) — > 7r+7r~?/;(2S). A fit to 
the mass spectrum with a single resonance yields a mass 
of (4324 ± 24) MeV/c 2 and a width of (172 ± 33) MeV, 
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FIG. 5: The measured CM energy dependence of the cross 
section (points with error bars) for e + e~ — > tv + tt~ ip(2S) af- 
ter background subtraction. The solid histogram shows the 
energy-dependent selection efficiency. 



TABLE I: Summary of main systematic uncertainties for the 
e + e~ — > 7r + 7r _ i/)(2S) cross section measurements. 



Source 


Systematic error 


Model-dependent acceptance 


±9.0% 


Tracking efficiency 


±7.6% 


B(ip(2S) -> 7T+7T- j/tp) ■ b(j/%i> -> e+e~) 


±3.5% 


Total 


±12.3% 



where the errors are statistical only. The structure in 
Fig. [3] has a mass that differs somewhat from that re- 
ported for the F(4260) in Ref. 0. However, the possi- 
bility that it represents evidence for a new decay mode 
for the U(4260) cannot be ruled out at this time. 
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